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ABSTRACT 


We demonstrate for a passive optical cavity containing a dispersive atomic medium, the 
increase in scale factor near the critical anomalous dispersion is not cancelled by mode broadening 
or attenuation, resulting in an overall increase in the predicted quantum-noise-limited sensitivity. 
Enhancements of over two orders of magnitude are measured in the scale factor, which translates 
to greater than an order-of-magnitude enhancement in the predicted quantum-noise-limited 
measurement precision, by temperature tuning a low-pressure vapor of non-interacting atoms in a 
low-finesse cavity close to the critical anomalous dispersion condition. The predicted enhancement 


in sensitivity is confirmed through Monte-Carlo numerical simulations. 
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1. INTRODUCTION 


The suggestion that the sensitivity of a passive optical cavity to a change in optical path length 
can be increased by placing an absorbing medium inside the cavity is counterintuitive. Near their 
resonance frequencies, not only do the atoms absorb energy decreasing signal-to-noise, the cavity 
line-width also broadens as a result of the associated anomalous dispersion, i.e., the white-light 
cavity effect. The combined effect would presumably reduce measurement sensitivity. However, 
this does not take into account the increased shift in the cavity mode frequency in response to a 
change in cavity length (or other external stimulus), i.e., scale-factor enhancement, which also 
occurs in the vicinity of such an absorption resonance. Recent publications have demonstrated that 
the scale factor is enhanced provided the dispersion is anomalous. Furthermore, the scale factor 
increases faster than the cavity line-width as the absorption and corresponding anomalous 
dispersion are increased, provided the cavity is under-coupled, resulting in a net increase in the 
scale-factor-to-mode-width ratio [1-3]. These increases can be particularly dramatic as the cavity 
approaches a critical anomalous dispersion (CAD) condition. Nevertheless, the increase in the 
scale-factor-to-mode-width ratio is accompanied by a substantial decrease in signal-to-noise as a 
result of medium absorption [4]. The question of whether the overall measurement sensitivity of a 
passive cavity can be increased has, therefore, persisted [5]. 

We address this question by simultaneously measuring the change in scale factor, mode width, 
and signal-to-noise as a low-finesse cavity is scanned in optical path-length and one of its modes 


passes through the D2 Fz = 2 to Fe Doppler-broadened resonance of a low-pressure vapor of non- 


interacting °’Rb atoms. By temperature tuning the atomic vapor, a pole in scale factor is observed 
near the CAD condition. A semi-empirical absorption/dispersion model [6], valid at low 
intensities, is used to extract the temperatures and scale-factor enhancements, while the change in 
mode-width and signal-to-noise are directly measured. The change in frequency resolution for the 
dispersive cavity in comparison with the corresponding empty cavity can then be inferred from 
these measured quantities. To check our assumptions and compute the frequency uncertainties 
directly, rather than through inference, we perform Monte-Carlo simulations using the cavity 
spectrum determined from the theoretical model as a master probability distribution function, from 
which sample photon distributions are randomly drawn. Finally, we determine whether an 
enhancement in sensitivity to changes in optical path length can occur under ideal (quantum-noise- 
limited) circumstances, 1.e., whether the increase in scale factor can be larger than the associated 
increase in frequency uncertainty. 

The low finesse of the cavity is particularly convenient in the experiments and simulations 
described above because it ensures: (i) the intracavity intensity is low enough for application of 
the absorption model, (ii) the cavity is not so under-coupled that the attenuation of the modes 
causes them to disappear entirely when they are tuned across the absorption resonance, and (iii) 
the width of the mode is sufficiently large that the Doppler-broadened absorption profile is not 
effectively uniform across it, i.e., the dispersion cannot be approximated as linear over the finite 
width of the mode. Whereas enhancement of the scale factor only requires linear dispersion, 
enhancement in the scale-factor-to-mode-width ratio requires higher-order dispersion. The latter 
results in mode reshaping which limits mode broadening and further enhances the scale factor, 
resulting in an increase in the scale-factor-to-mode-width ratio. Without this mode reshaping, the 
scale factor and mode width would increase by the same amount and their ratio would not increase 


above unity. Mode reshaping is, therefore, crucial to achieving an increase in sensitivity. The 


relevant parameter to obtain mode reshaping is the ratio of the cavity-mode-width to absorption 
resonance width, which should not be too much smaller than unity. 

While the low cavity finesse certainly makes our experiments easier, it is by no means 
fundamentally important to obtaining an enhancement. Rather, the crucial requirement is that the 
dispersive resonance width and cavity linewidth be comparable. For the relatively wide Doppler- 
broadened resonances used here, this condition happens to correspond to a cavity with a low 
finesse. However, if the atomic resonance used for the dispersion were narrower, the required 
finesse would be correspondingly higher, i.e., high-finesse cavities can also be enhanced but 
require narrower absorption features (so that the dispersion cannot be approximated as linear). 
Therefore, our conclusions apply equally well to the case of high-finesse cavities, which are more 
technically relevant for devices owing to their higher absolute sensitivities. 

Previously, we measured the cavity scale-factor enhancement as the temperature of the 
intracavity dispersive atomic medium was tuned, and demonstrated how temperature fluctuations 
in the atomic medium limit the scale-factor stability [6] and, therefore, the measurement precision. 
In the present work, we limit our discussion to the quantum noise limit (QNL) by ignoring classical 
noise, and therefore neglect the effect of these temperature fluctuations as they are not relevant in 
the QNL. Furthermore, to address the question at hand, we extend the previously-obtained 
experimental results, which were limited to scale factor, by providing additional measurements of 
the cavity mode width and signal-to-noise while temperature tuning. 

An enhancement in optical cavity sensitivity using anomalous dispersion could have 
applications ranging from detection of gravity waves [7-10], tests of general relativity [11], optical 
communications schemes [12-14], enhanced strain and displacement sensing [15], wideband 
coherent perfect absorbers [3, 16], and increasing the precision of optical gyroscopes [1-3, 6, 11, 


17-20]. Moreover, the fast-light enhancement is a particular example of a broader phenomenon 


that can be found in any physical system that possesses an exceptional point [3]. Such points are 
commonly found in non-Hermitian systems such as coupled oscillators having different loss rates. 
The critical anomalous dispersion is, in fact, an exceptional point, arising from the coupling of 
atomic and cavity resonant modes. Therefore, the enhancement in sensitivity described herein 
could ultimately benefit applications that extend far beyond those that rely on optical-cavity-based 
sensing. 

In section 2 of the text, the enhancement in the quantum-noise-limited (QNL) precision of a 
dispersive cavity is derived in terms of the scale factor, mode width, and signal-to-noise. In section 
3 the scale-factor enhancement in transmission and in reflection is derived for a cavity containing 
an atomic medium, and in section 4 an applicable model [6] for the absorption and phase shift of 
the medium is presented. Sections 5 and 6 discuss the experimental results, while section 7 


describes the results of the Monte-Carlo simulations. 


2. QUANTUM-NOISE-LIMITED PRECISION OF A DISPERSIVE CAVITY 


A useful metric for a sensor utilizing the measurement of cavity mode frequency displacement 
for a given change in cavity optical path length, e.g., in an optical ring cavity gyroscope, is given 


by 


Y= q a (1) 


where w, is the peak frequency of the gq" cavity mode, w|, is its peak frequency at the initial or 
reference cavity length, and éq@ is the measurement uncertainty in the peak frequency. Thus, if 
mode peak displacement w, — a! can be increased without a concomitant increase in dw, the 
sensitivity of the sensor is enhanced. The minimum resolvable input to the cavity may then be 


found from Eq. (1) by setting vy =1. 


It follows that the sensitivity enhancement (or reduction) of a dispersive cavity, compared to 


the corresponding empty cavity, can be written as 


—a) (e) 
pa = fe (2) 
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We denote the mode frequency, initial mode frequency, measurement uncertainty, full-width-at- 


r(e) 
’ 


half-maximum (FWHM) linewidth, and signal-to-noise ratio of the empty cavity by a” » W, 


© and SNR, respectively. The corresponding quantities for the dispersive cavity will 


6a, y 
not be superscripted by (e). The cavity scale-factor enhancement can then be defined as 
S(w,)=dw,/dw,’. Hence, for a dispersion-enhanced cavity operating in the vicinity of an 


absorption resonance, a small change in the cavity’s optical path-length will result in a mode peak 


displacement 

ri (e) r(e) 
w -o), =S(w —W ). (3) 
Therefore, regardless of the source of the uncertainties, the sensitivity enhancement can be written 


ba 
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(4) 


An equally valid way to arrive at Eq. (4) is to set y=1 and realize that an error, dw, in a 
dispersive-cavity is equivalent to éw/S in the corresponding empty-cavity. Note that the 
reduction in the minimum measurable cavity mode frequency is simply given by 1/€. 

Now let us assume that all classical sources of error in the cavity can be eliminated, e.g., due 
to temperature and mechanical fluctuations, leaving only quantum-mechanical photon shot noise. 
It has been reported [21, 22] that the frequency error for an empty (unidirectional) passive cavity 
in this QNL is 


(e) yo 
6@ont 7 SNR . (5) 


If we assume Eq. (5) also applies to a cavity containing a dispersive medium, then édw,, * y/ SNR 


. Accordingly, Eq. (4) becomes 


ONL : 6 
¢ Wr (6) 


where M =SNR/ SNR is the normalized signal-to-noise and W=y/y is the normalized 
cavity linewidth. Hence, if w e could be measured to the same QNL precision as a”, the 
enhancement would simply be ¢°” = S . However, because of mode broadening and attenuation, 
6Moy, 18 typically larger than 60 our . As we shall see, however, it is not larger by a factor of S. 
Mode broadening and attenuation reduce, but do not cancel out, the increase in scale factor, 
resulting in an overall increase in the QNL sensitivity. Note that if, in Eq. (6), we take 
SNR = SNR , i.e., M =1, then we obtain the sensitivity enhancement factor that we used in our 
earlier work [2], namely £°” =S/W. Also, note that while each of the quantities in Eq. (6) are 
frequency dependent, we will primarily be concerned with their values at a particular frequency 
w,, defined as the frequency where the scale factor enhancement is maximized (near the 
absorption peak). The QNL enhancement, scale-factor enhancement, normalized signal-to-noise, 


and normalized mode width at w, will be denoted by £°", S 


»> M,, and W,, respectively. 


In the experiment outlined in section 5, we use Eq. (6) to evaluate the peak QNL sensitivity 
enhancement, £°", by measurement of S,, M,, and W,. The advantage of this approach is that 
one can predict the QNL enhancement without having to directly measure frequency errors, 1.e., 
perform a full noise analysis. However, a critical assumption behind Eq. (6) is that Eq. (5) should 
apply to the dispersive cavity in the same proportion that it applies to the empty cavity. This 
assumption neglects the effect on the uncertainty of the change in mode shape that occurs as the 


CAD condition is approached (the mode becomes more flat-topped). In section 7, on the other 


hand, the assumptions leading to Eqs. (5) and (6), are abandoned, and we take Eq. (4) as our 


starting point for finding the QNL enhancement. The frequency errors éw,y,, are calculated 
numerically via Monte-Carlo simulations using the cavity transmission spectrum, see Eq. (10), as 
a master photon probability distribution function, from which sample photon distributions are 
drawn. Because they are drawn randomly, the statistics are Poissonian, and therefore the computed 
frequency errors are the QNL quantities of interest. The importance of this numerical procedure is 
that it bypasses the assumptions behind Eq. (6), providing a more direct indication of the QNL 
enhancement that takes full account of the change in mode shape near the CAD condition. We find 
that the RHS of Eq. (6) must be modified by a constant factor K that depends on the particular 
form of averaging employed in finding the distribution peak. 

We assume that the source of the quantum noise in our experiments is detector shot noise and 
any additional sources of quantum noise, e.g., due to spontaneous emission in the medium, are 
negligible. This is a good assumption for passive cavities because medium absorption and optical 
pumping are much weaker, and therefore the relative amount of spontaneous emission to laser 
light is much smaller, than for the case of a laser (where spontaneous emission noise cannot be 
neglected). The signal-to-noise for shot-noise-limited direct detection of photons on a photodiode 


is 


SHR be IN. (7) 
S 


where JN, is the number of signal photons per second arriving at the detector. In this case the 
signal is associated with the resonance of a cavity. Eq. (7) can be generalized to account for cavity 
fringe baselines by substituting N,—N, for N, in the numerator and N, +N, for Ny (variances 
add) in the denominator of Eq. (7), where N, is the number of photons per second arriving at the 


detector for the baseline (away from the cavity resonance). This is equivalent to substituting 


(N, —N,) /(N, +N,) for N, on the RHS of Eq. (7). The signal-to-noise can, therefore, be 


obtained from transmission or reflection spectra by using the relation 


(Ns—-N,) _( m_\P, _( 2V’)P, (8) 
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where m=|N, —N,|/N, is the fringe modulation depth, V is the fringe visibility, and P, is the 
baseline power, and the positive (negative) sign is used in transmission (reflection). Note that the 
factor P,/ha is the same for the dispersive and empty cavities and hence factors out in the 


determination of the ratio, MM . Therefore, 


oe se) (e) 
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where the relation m=2V/(1+V) has been used to transform to fringe visibility. Again, the 


superscript (e) is used to denote the empty-cavity fringe modulation depth and visibility. 


3. SCALE FACTOR OF A DISPERSIVE CAVITY 


The complex reflection and transmission coefficients for a Fabry-Pérot or two-port ring cavity 
containing a dispersive atomic medium at temperature T can be written as 


1 i _ 2(a, Tye 
I, J- 2(a, Te % 


p(,T)= (10) 


and 


~ ior ,/2 
#(o,T) = hts V8(o,T)e See an) 
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respectively, where 7, is the round-trip time of the empty cavity. The net electric-field gain per 
round trip is g(w,T)=7n74,7,,(@,T), where r,, and f,, are the real-valued reflection and 


transmission coefficients of the mirrors, respectively, a, accounts for other frequency-independent 


round-trip losses in the cavity, and 7,,(w,T) =7,,(w,T)exp[iP(w,T)] is the complex round-trip 
transmission coefficient of the atomic medium. 

The resonance condition is then determined by finding the frequencies w, where the derivative 
of the transmission lH, 7) or reflection |p(@,7)| with respect to w goes to zero, i.e., from the 
transcendental equation 


w,T,.+P(w,,T)+ F (w,,T) =27q = wr, , (12) 


where q is the cavity mode number, ®(w,7) is the effective round-trip phase-shift of the medium, 
and F(w,T) is an additional phase factor that arises from the reshaping of the mode by the medium 
absorption. The scale-factor enhancement is then obtained by taking the derivative of the 
dispersive cavity mode frequencies with respect to the empty cavity mode frequencies, 


$(w, ,T) = : 
w,,T)= =— 
do? A, (w,,T)+T. 


qd cav 


en’ (13) 


where 71,(w,T)=1+(1/7,)d B/ dw is the effective group index and T,,,(@,T) =(1/7,)d F/dw 
is the additional dimensionless time delay associated with absorption. Note that the scale factor 
has a pole at a critical anomalous dispersion, and is enhanced when the sum of the group index 
and this additional cavity delay time is less than unity, 1.e., the sum of the slopes or net time delay 
contributed by these two terms must be negative for an enhancement to occur. 


The functional form of F in transmission differs from that in reflection, which leads to different 


scale factors for transmission and reflection. In transmission 


(14) 


F(@,T)=—sin™ ae ake . 


2g(w,T) f,(w,T)r, 


where g'(w,T) =dg / da. On the other hand, in reflection 


Fooit)=—sn"| A@.T) Jose Co) } (15) 


{A?(0,T) + B?(o,T) JA2(w,T) + B2(w,T) 
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where 
A(w,T) = g'(w.T)(7? + (wT) ) 
B(w,T) = g(o,T)A,(@,T)T, (77 — g(@,T) ) (16) 
C(w,T) = g(w,T)g(@.T)(77 +1). 
Note that while the phase shift ®(w,,7) influences the mode frequencies by mode pushing, 
F (@,,T) does so instead by mode reshaping. When the cavity mode is narrow in comparison with 
the atomic resonance, 1.e., in the linear-dispersion approximation, the absorption variation across 
the mode and resultant mode reshaping can be neglected, and we obtain S(T) =1/n 7) [11, 20]. 
The cavity linewidth iS given, within this approximation, by 
y(T) =(4/ A, (T)r,)sin | Jd=8)(2)/20+ g,(T)) . Therefore, the normalized cavity linewidth 


can be written as 


(T)_ wu) 
W (T) = oo - _ (T) (17) 


where u(7) takes into account the change in the resonant value of the net round-trip field gain, g, 
with temperature. The function u(T) approaches unity at low temperatures and exceeds unity at 
all other temperatures. Therefore, W(T) is always greater than 1/n,(T). Hence, if the linear- 
dispersion approximation were valid for all scale factors, there could be no enhancement of the 
scale-factor-to-mode-width ratio or QNL sensitivity, because broadening of the mode would 
always cancel the increase in scale factor. This cannot happen at all scale factors, however, because 
the broadening of the mode ensures that at some point the linear-dispersion approximation must 
be violated, in particular as the width of the mode approaches that of the medium resonance. 
Therefore, the scale factor and critical anomalous dispersion are more properly determined by Eq. 
(13). In addition, the mode width does not increase as fast as predicted by Eq. (17), but is instead 


clamped by the finite medium absorption width. The increased scale factor and decreased mode 


11 


width, in comparison with the linear-dispersion prediction, work together to increase the scale- 
factor-to-mode-width ratio and QNL sensitivity. 

According to Eq. (15) there is another way to obtain F(w yl) =0, namely when 
A(w,T) =C(@,T). The ratio A(w,T)/C(@,T) thus determines whether mode reshaping from 
F (@,,T) will add to or subtract from the mode pushing due to ®(w,,7), i.e. whether the slope 
T.,,(@,T) will be negative or positive, respectively. Enhancement of the scale-factor-to-mode- 
width ratio (and, therefore, of the QNL sensitivity) requires that the mode reshaping adds to the 
mode pushing, ie., T..(@,T)<O or A(w,T)>C(@,T). In reflection this only occurs for 
g(w,T) <r, or 74,7,,(@,T) <4, i.e., when the cavity is under-coupled. This requirement is 
always fulfilled for symmetric cavities (7,=7,). Single-port (all-pass) cavities like the one 
presented in section 5, on the other hand, require a,T,,(@,7) <7,. For over-coupled cavities the 
scale factor increase is more than offset by increased mode width, resulting in a decrease in the 
scale-factor-to-mode-width ratio. In transmission, on the other hand, only one term appears in Eq. 
(14), and the mode reshaping always adds to the mode pushing. 

Note Eq. (10) assumes that any mirror losses, other than reflection, in the transmission through 
the input mirror are negligible, i.e., 7° +t) =1. We can take any additional loss into account by 
assuming the transmission coefficient of the first mirror is t7 =1—77 —b; , where b, accounts for 
this additional loss. The form of Eq. (11) and subsequent equations for transmission are unchanged 
by the added mirror loss. The equations for reflection, on the other hand, require substituting 
g(1—b;) for g in the numerator of Eq. (10) and wherever it, or its derivative, appears in 


Eq. (16). 
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4. ABSORPTION AND DISPERSION MODEL 


For a given cavity, the functions B(w,T) and F (w,7T) are uniquely determined by the medium 
absorption coefficient a(w,T). As we discussed previously in detail [6], a semi-empirical model 
that takes into account all the Zeeman-degenerate Doppler-broadened hyperfine transitions of the 
87Rb D> transition can be used to model the absorption provided the intensity is sufficiently weak 
that optical pumping, saturation, and power broadening are all negligible. Under these conditions 
the atoms remain in thermal equilibrium with level populations determined by the Boltzmann 
distribution. The absorption coefficient is then simply a weighted superposition of the various 


hyperfine transitions, 


(18) 


. 


u ; WW, ; a 
a(w,T) ouDZ0.jo | r (T) 


where w, , and s(2, j) represent the frequencies and strengths of the various hyperfine transitions 
from the F,=2 ground state, respectively. The sum is over the excited states Fe = j. The 
temperature dependence of a(w,7T) is determined by that of the 1/e Doppler width 
r,(7)= k[2k,T /m and envelope function a,(T)=C-N(T)/1T',(T), where C is constant to a 
good approximation. We assume the number density N(7) follows the empirical relation 
N(T) = (133.3/k,T) 1008" | where A=4.312, B=-4040, and T is specified in Kelvins 
[23]. Note, the coefficient a,(T) represents the peak absorption coefficient only when the 
absorption can be properly described by a single Gaussian. Assuming linear polarization such that 
only transitions between states having the same magnetic quantum number are allowed, and 
ignoring Zeeman pumping so that the distribution of ground-state magnetic sublevels is uniform, 
the coefficients s(2, j) for the allowed transitions are s(2,1) =0.05, s(2,2)=0.25, s(2,3)=0.7, 


and C =(5/8)o07 ViT/2=4.1x10°m?/s, where T is spontaneous decay rate from the excited 
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state and 07 =1.938x10m” is the total on-resonant absorption cross-section for 2-polarized 
light. 

The medium effective phase-shift ®(w,7) can then be determined from a Kramers-Kroénig 
(K-K) relation by taking the Hilbert transform of the medium transmission coefficient 
T,,(w,1) =exp[—a(w,T)C/2], where ¢ is the length of the atomic medium. For a Gaussian 
absorption profile this transform has the simple analytic form 
H {exp[—aexp (—x° / b*)]} = ia-erf (ix/ b)-exp(—x° / b*) . The phase shift is found by applying this 


transformation to each of the individual hyperfine transitions, i.e., 


(0,7) = POLY 502, jen (5-2 on (i ) (19) 


P77) P,() 


j=0 


5. DESCRIPTION OF EXPERIMENT 


The experiment is the same as that we reported previously [6], except here we also measure 
mode widths and signal-to-noise, and our approach to the analysis of the results differs (see section 
6). To summarize, a linearly-polarized external-cavity diode laser at a wavelength of 780 nm was 
used to scan over the modes of a 40 cm ring cavity, containing a 2.5 cm isotopically-enriched ®’Rb 
quartz vapor-cell enclosed in a temperature-stabilized aluminum oven [6, 24]. An all-pass cavity 


configuration (7, = 1, |T(@, T)| =0) consisting of three high reflectance mirrors and a 90:10 cube 


beam-splitter for the input/output coupler was used, as shown in Fig. 1. The center frequency of 
the laser was adjusted to coincide with the Doppler-broadened F, =2— F, resonance and the 
laser was scanned over several cavity free spectral ranges. A Michelson interferometer with 
unequal arm lengths and a second *’Rb cell in a counter-propagating pump-probe saturated- 


absorption spectrometer were used as frequency references, and a reference detector was used to 
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factor out variations in the laser intensity during the spectral scan. The incident beam was s- 
polarized, and the detuning between the cavity and atomic resonance was varied by an intracavity 


liquid-crystal variable retarder whose slow (tuning) axis was aligned with the polarization. 


Fig. 1 (Color online) Experimental Setup. BS = Beam-Splitter, ISO = Isolator, LC = Liquid Crystal Variable 
Retarder, IC = Input Coupler, SA = Saturated Absorption Spectrometer, Ref = Reference Detector, MI = 
Michelson Interferometer, R = Cavity Reflection. 


The temperature of the atomic vapor cell was varied from 36°C to 42.5°C through the critical 
anomalous dispersion (which occurred near 42°C). At each temperature, reflection spectra were 
recorded at a variety of detunings as a selected cavity mode was tuned across the atomic resonance 
via the variable retarder. Representative experimental spectra are shown in Fig. 2. A larger number 
of spectra were recorded close to the resonance by applying a nonlinear voltage step to the retarder. 
An automated peak finding program was used to obtain the frequencies, widths, and modulation 
depths of modes near to and far away from the resonance, resulting in plots of the mode detunings 
for the dispersive cavity vs. those for the empty cavity as shown in Fig. 4(a). In this figure, each 
data point represents a different liquid crystal voltage and the detunings of the dispersive and 


empty cavities are converted into linear FSR units, ie., the detunings are defined as 
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A=(w,-@,)t,/2m and &=(w\° —w,)t,/27, respectively, where w, is the frequency of the 
scale factor maximum, which, owing to the asymmetry of the absorption profile, is shifted slightly 
from the absorption peak. The subscript on the quantities S,, M,, W,, and f° refer to @,, i-e., 


Sy =S(@, = @p)), ete. 


6. EXPERIMENTAL RESULTS 


In this section the relevant parameters S,, M,, and W, are obtained from the experimental 


O2 


reflection spectra for evaluation of £2” via Eq. (6). Representative spectra are shown in Fig. 2 for 


the case where the cavity is held near the CAD condition. Several points are worth remarking 
upon: (7) the effect of detuning the mode is greater for the mode nearest to the atomic resonance ( 
A>6), 1e., the scale factor is enhanced for this mode, (ii) this mode nearest resonance also 
broadens and is attenuated, (iii) the mode broadening is limited by the width of the atomic 
resonance, and (iv) the mode amplitude is never fully extinguished even at the CAD condition. 
Therefore, even without explicit evaluation of Eq. (6), simple observation of the cavity spectra 
suggests that a QNL sensitivity enhancement is possible, because at the CAD condition the scale 
factor enhancement approaches infinity, while the mode width and amplitude always remain finite 


and non-zero. 
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Reflectance, R 


=D -1 0 1 
Laser Detuning (units of FSR) 


Fig. 2 (Color online) Representative reflection spectra data taken near the critical temperature (sensor 
temperature 314.25K, fit temperature 314.77K), at two different values of the empty-cavity detuning 6 (solid 
curve is for 6 =0 ). A saturated absorption spectrum (top) is also shown for reference. As the mode approaches, 
it is pushed away and reshaped by the atomic resonance, resulting in a larger dispersive-cavity detuning A. 
Consequently, the scale factor increases. The mode also broadens and is attenuated, but note the broadening is 
limited to the Doppler-broadened resonance width. As a result the attenuation is also less than expected. The 
on-resonance normalized mode width and signal-to-noise, W, and M,, respectively, are calculated from 
measurements of the FWHM mode width and modulation depth (the dip depth divided by the fringe baseline), 
taken both on-resonance and away from the resonance, as shown. 


The point of the current work is not to validate the theoretical model for S, vs. T (as it was in 
our previous work [6]), but rather to use it to extract the maximum scale-factor enhancements by 
fitting to the experimental data. The following procedure was used to determine the scale-factor 
enhancements: For each spectrum, the FSR was measured. The average value of 1000 such 
measurements was 703.0 MHz with a standard deviation of 3.1 MHz, i.e., uncertainty of 0.1 MHz. 
The value of 6 could not be directly measured, but was inferred by using a mode located 
sufficiently far from the influence of the atomic resonance and adding the appropriate number of 
FSR to its frequency. At each temperature, the two parameters, 7, and a,, were found by a 
nonlinear least-squares fit of the theory to the part of the reflection spectrum that was away from 
the atomic resonance. When input mirror losses (other than reflection) were ignored, the average 
values of these parameters were a, = 0.674 and 7, =0.946 with standard deviations of 0.044 and 


0.006, respectively (over 22 measurements). When these losses were considered, the average best 
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fit values were a, =0.672, 7, =0.917 and b, =0.209, with standard deviations of 0.020, 0.040, 
and 0.094, respectively (over 6 measurements). For comparison, direct measurement of the mirror 
transmission and reflection coefficients yielded 7,=0.903 and b, =0.260. While the better 
agreement suggests that mirror losses should be considered, the inclusion of these losses ultimately 
makes no difference in the determination of S,, M,, or W,, because the ratio of the cavity-mode- 
width to absorption resonance width is unchanged. From these parameter values it is obvious that 
the cavity was under-coupled, as required to obtain an enhancement in the scale-factor-to-mode- 
width ratio. 

The largest value of the input intensity was 4.2 mW/cm/? for these measurements (the typical 
value was 1.5 mW/cm?). However, it is the steady-state intracavity intensity that is important for 
determining the applicability of the absorption model. The average cavity finesse for 10 different 
temperature measurements was 7.6 with a standard deviation of 1.3. The intracavity intensity at 
the peak of a mode on-resonance with the atoms was always less than 20 pW/cm’, far below the 
lowest saturation intensity for the hyperfine transitions with z-polarized light (3.6 mW/cm’ for the 
| 2 to F, =3 transition [23]). The low finesse and large on-resonance absorption self-limit the 
steady-state circulating intracavity intensity, enabling application of the model. 

Each set of data in Fig. 4(a) consisted of 500 data points, obtained from spectra of the sort 
shown in Fig. 2. From each data set the scale-factor enhancement was determined by performing 
a nonlinear least-squares fit to the theory presented in sections 3 and 4 (Eq. (12) was the fit 
equation) with temperature as a fit parameter. The fitted temperatures are compared with those 
measured by a platinum RTD sensor mounted on the outside of the aluminum oven in Fig. 3. The 
strong correlation appears to confirm the theoretical model is valid for the intensities used in the 
experiment. However, systematic errors are likely present that weaken this argument. The most 


important of these are: (i) the sensor temperature may not accurately represent the temperature of 
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the vapor (shifts the abscissa values), and (i7) the values of the A and B coefficients (see section 4) 
in the empirical relation for the number density may need to be adjusted (shifts the ordinate values), 
rather than relying on the published values [23] for an infinite liquid reservoir as we do here. We 
have observed that these coefficients vary over time for a given vapor cell. Estimated upper bounds 
for these systematic errors as determined from ancillary measurements [6] are +2K for the abscissa 
and +2.5K for the ordinate. What can be concluded, though, is that the optical detection/fitting 
procedure provides a more precise (if not necessarily accurate) temperature measurement than our 
temperature controller/sensor system. The estimated single-measurement temperature uncertainty 
of the RTD sensor due to temperature fluctuations and limited readout resolution was 150 mK 
(abscissa) (this value can be decreased through additional passive controls and/or by making 
multiple measurements), an order of magnitude larger than that obtained from the fitting 


procedure, 15 mK (ordinate). 


Model Temperature (K) 


310 312 314 316 
Sensor Temperature (K) 


Fig. 3. (Color online) Temperature determined from the model vs. temperature measured by an oven-mounted 
sensor. The dashed line is a linear fit, weighted by the errors in both coordinates. The slope and intercept of 
this line were 1.004 and -0.926 K, respectively, with a correlation coefficient of 0.917. 


The temperatures resulting from the fitting procedure and maximum on-resonance scale factors 


are plotted in Fig. 4(b). Note that a pole occurs in the scale factor corresponding to a critical 
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anomalous dispersion of 7, = 0.483 and critical temperature of T, =314.82 K. This is far from 
the values predicted in the linear-dispersion approximation (7, =0, T, =322.13K ). The highest 
scale-factor enhancement recorded was S, = 363 (with errors of +1175 and —157). Note the error 
increases quadratically with the scale factor, and at some point will even exceed the value of the 
scale-factor enhancement itself. The error bars in Fig. 4(b) represent the uncertainties derived from 
the fit parameters assuming a 5 MHz error in the determination of A and a 6 MHz error in 6. 
These errors primarily result from cavity noise during the spectral scan. (The FSR measurement 


makes the error in 6 slightly larger than that of A .) 
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Fig. 4. (a) Scale factor plots near and far from the critical temperature (314.82 K). The detunings, A and 6, 
are in units of FSR. The temperatures are obtained from the nonlinear-least-squares fits (solid curves). The 
curves have been shifted slightly for clarity. (b) Scale-factor enhancement vs. best fit temperature. The points 
(i) and (ii) represent the maximum scale factor slopes obtained from the corresponding data sets in (a). The 
dotted curve shows the expected results from the model when average values are used for 71 and a1. 


From the cavity spectra, the on-resonance mode width and signal-to-noise were also measured 
and plotted normalized to their empty cavity (off-resonance) values in Fig. 5 (the data points). The 
temperatures of the data points in Fig. 5(a) are those from Fig. 4(b), i.e., they were computed using 
the full theory. In Fig. 5(b), on the other hand, the temperatures were computed by using the linear- 


dispersion approximation to find the temperatures that correspond to the scale-factor 
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enhancements in Fig. 4(b), leading to an average shift of +6.6K, slightly less than the shift in the 
critical temperature (+7.3K). The theoretical curves were determined numerically from spectra 
generated from the full theory (Fig. 5(a)), as well as analytically using the linear-dispersion 
approximation (Fig. 5(b)), i.e., Eq. (17), using the average values of 7, and a,. Note that for both 
theoretical curves, M, and W, approach unity at low temperatures, as expected when the medium 
does not absorb. At higher temperatures, however, the two theories diverge. In the linear- 


dispersion approximation, W, shares a pole with S,,i.e., at 1/n, , whereas in the full theory W, 


0°? 
does not increase as fast and remains finite, because the finite width of the absorption profile 
clamps the increase in mode width. Moreover, M, converges to a finite value, whereas it vanishes 
in the linear-dispersion approximation. Note that the linear-dispersion approximation does not 
match the data in Fig. 5 (b). The deviation occurs because, as discussed in section 3, the linear- 
dispersion approximation is not valid because the mode width is not sufficiently narrow compared 
to the absorption resonance (see Fig. 2). As a result, the full theory matches the data better for both 
W, and M, in Fig. 5 (b). However, the measured value of M, is still consistently higher than the 
full theoretical prediction at all temperatures (by a constant shift of 0.078). This offset could be 


the result of a small systematic error in the determination of the modulation depth for the off- 


resonant mode. 
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Fig. 5. (Color online) Factors that reduce sensitivity: on-resonance normalized mode width (top) and signal-to- 
noise (bottom) vs. vapor-cell temperature. For the data points, the temperatures were determined from the 
nonlinear least squares fits, whereas W, and M, were determined by direct measurement from cavity spectra. 
The data points are compared with the (a) full theory and (b) linear-dispersion approximation, using average 
values for ri and ai. The temperature of the pole in W is indicated by the vertical line in (). 


In Fig. 6(a) the QNL sensitivity enhancement obtained from Eq. (6) is plotted against the 
temperatures obtained from the model. As expected the enhancement is less than unity below the 
critical temperature, but then rises above unity closer to the critical temperature. Again, the data is 
consistently higher than the theoretical prediction owing to the offset in M,. The largest value 
measured was 18.8 (with errors of +61 and —8.5). Note at low temperatures the error is roughly 
constant as it is determined primarily by the measurement of W, and M,, but as the critical 
temperature is approached the error increases dramatically due to the increasing contribution from 
the scale factor measurement. 

A comparison between the full theory and linear-dispersion approximation is shown in Fig. 
6(b). Note in the linear-dispersion approximation there is no enhancement in QNL sensitivity, 
because S, and W, both increase as 1/n, . Therefore, the scale-factor pole at 322.13K does not 
appear in the overall QNL enhancement. In the full theory, on the other hand, the pole in scale 
factor persists in the QNL sensitivity. Again, there are three reasons for this: (7) the scale factor 


increases faster than 1/n, due to mode reshaping, (ii) the increase in mode width is limited by the 
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finite atomic absorption width / higher-order dispersion, and (iii) the decrease in signal-to-noise is 


limited because the cavity is not as far from critical coupling as it would be in the linear-dispersion 


approximation. 
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Fig. 6. (Color online) Predicted QNL sensitivity enhancement vs. vapor cell temperature. (a) Experimental 
results demonstrating that the enhancement rises above unity as the critical anomalous dispersion is 
approached. The solid curve represents the theoretical prediction. (b) Comparison of the full theory (solid 
curve) and linear-dispersion approximation (dashed curve). For each case, an arrow indicates the temperature 
where the scale-factor pole occurs (314.82K for the full theory and 322.13K in the linear approximation) 
corresponding to the critical anomalous dispersion. The pole does not show up in the QNL enhancement in the 


linear-dispersion approximation, but does appear in the full theory. 


7. MONTE-CARLO SIMULATIONS 


Recall that a crucial assumption in the derivation of Eq. (6), and therefore in the calculation of 
the QNL enhancement in the previous section, was that the effect of the change in mode shape on 
the frequency uncertainty of the mode could be neglected, even as the CAD condition is 
approached. We clearly see in Fig. 2, however, the modes becoming more flat-topped as they 
approach the atomic resonance. Therefore, in this section we drop the assumption behind Eq. (6), 


and compute the frequency errors for the dispersive and empty cavities through numerical Monte- 


Zo 


Carlo simulations. This procedure provides a more direct indication of the QNL enhancement 
through use of Eq. (4), which takes account the change in mode shape near the CAD condition. 

In this section it will be more straightforward to work with transmission, rather than reflection. 
We begin by assuming that the cavity transmission spectrum T(@) = |z(@)| represents a 
probability distribution for the occurrence of photons at frequency @, and then draw random 
sample distributions from this underlying master distribution. An alternative approach would be 
to start with a master transmission curve and simply add Poissonian noise to it to generate the 
sample distribution [25], but we found this procedure yielded similar results and so do not present 
it here. 

The central limit theorem states that the mean w of a sample probability distribution will be 
normally distributed regardless of its particular shape. As a consequence, the standard deviation 
of the mean will vary with the number of samples N according to o,,=0/ JN , where o is the 
standard deviation of the original sample distribution. For unimodal distributions that are not 
skewed, the mean provides an estimate of the peak of the underlying distribution. It follows that, 
in the quantum noise limit, the frequency error associated with the determination of the 


transmission peak is 
6@on, =a: (20) 


where SNR = JN is the signal-to-noise ratio at the detector, N is the number of photons, o is the 
standard deviation of 7(@) , and we have assumed for convenience that detector efficiency is unity. 
A requirement for Eq. (20) to be valid is that the SNR include all the photons in the sample 
distribution. This definition for the SNR is slightly different than in the discussion leading up to 
Eq. (9) in section 2, which is used to calculate M from our experiment, where only photons at the 


frequency peak were considered to contribute to the signal. We point out that the frequency range 
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of the sample photon distribution depends on the particular experimental design (the frequencies 
of passive gyros, for example, are locked to the cavity resonances), and this distinction is negligible 
when the range is small. Moreover, while the absolute value of the SNR obtained from the two 
definitions can be quite different, the normalized quantity M is not significantly modified and 
varies similarly with cavity scale factor (in the limit of high finesse the dependency of M on scale 
factor is the same), irrespective of the choice of definition. 

We have verified the validity of Eq. (20) for empty-cavity mode spectra over a large range of 
SNR by performing Monte-Carlo simulations to calculate 0, for a large number of random photon 
distributions. A practical difficulty arises because these spectra are not unimodal, and do not 
possess a standard deviation that converges, at infinite sample range. Spectra can be truncated to 
a single free spectral range such that Eq. (20) applies. Even so, the standard deviation will increase 
with sample range because cavity spectra are approximately Lorentzian at high finesse, and 
therefore always have significant probability density in their wings. It therefore makes sense in 
this case to use a range-independent parameter such as the full-width-at-half-maximum (FWHM) 


y to characterize the mode width, 1.e., 


Y 
60,,, =C——, 21 
ONL SNR (21) 


where the constant C=o/y is dependent on range and cavity finesse. Note that the approximate 
relation of Eq. (5), relied upon extensively in the gyroscope literature, has been replaced by an 
exact relationship, i.e., Eq. (21). Notably, for high-finesse cavities (such as gyroscopes) C can 
differ significantly from unity. 

An additional difficulty applies for cavities that contain a dispersive medium. In this case the 
proportionality constant in Eq. (21) also depends on the scale factor S, because the mode broadens 


and therefore takes up a larger portion of the data range. In this case, Eq. (6) can be rewritten as 


2a 


5 (e) . 
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cone = 
where K =C /C. Moreover, because the modes can be skewed when the cavity mode is shifted 
slightly from the medium resonance frequency, the distribution mean does not accurately represent 
the mode peak. Eqs. (20)-(22) are not, in general, valid under these circumstances. Nevertheless, 
the peak can still be found by curve fitting, other averaging techniques, or simply by locating the 
mode of the distribution. This approach has the advantage that, to the extent they may involve less 
averaging, these alternate peak finding methods can be less dependent on data range. The 
distribution mode, for example, is independent of data range for sufficiently large SNR. In fact, 
one would expect the mode method to produce a limiting worst-case-scenario result as it uses the 
least amount of averaging, whereas the mean method represents the best-case limiting behavior, 
ie., they represent cases of maximum and minimum amounts of averaging, respectively. Hence, 
the value of ¢°” that can be achieved in a given situation depends on the method used to identify 
the peak of the photon distribution, and how much averaging is involved in that method. 
Therefore, we performed Monte-Carlo simulations using both these limiting methods. For each 
data point, the error 6@,,, was obtained by taking the standard deviation of the mean or mode 
values obtained from 250 different photon distributions. The data range was limited to a single 
empty-cavity FSR, centered on the transmission peak. The same cavity parameters were used as 
in section 6, namely a, =0.674, 7, =0.946, and FSR of 703.0 MHz. In addition, it was assumed 
that a second identical mirror with r, =0.946 was used as an output coupler, so that cavity 
transmission could be monitored (rather than reflection as in section 6). The finesse in this case is 
6.5, slightly lower than that used in our experiments in section 6. Photon distributions were 
produced by mapping a uniform random distribution onto the cavity transmission function through 


inverse transform sampling, i.e., by taking the squared-magnitude, integrating, and inverting Eq. 
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(11). The maximum scale-factor enhancement was varied from S,=1 to S, =344 by changing 
temperature, and signal-to-noise was varied by changing the number of photons. In addition, as 
the scale factor increased, the photon number (and SNR) was reduced to account for attenuation of 
the mode as a result of medium absorption. 

In Fig. 7, the ratio K =C/C is plotted in transmission as the scale factor is varied. Recall 
that for the mean method, knowledge of C and K allows the frequency errors and QNL 
enhancement to be determined apriori by application of the central limit theorem (Eqs. (21) and 
(22)). Hence, for the top curve in Fig. 7, K is calculated directly from the values of o and y for 
the master photon distribution, and is therefore independent of SNR. Note also that the value of K 
asymptotes to a constant at high scale factor. This convergence of K occurs because, as the scale 
factor increases, the broadening of the mode is eventually clamped by the finite absorption 
resonance width, i.e., by higher-order dispersion. At higher values of the finesse, K converges to 
a larger value. For the bottom curves, on the other hand, the value of K was determined 
aposteriori, i.e., by using Eq. (22) to back out the value of K after computation of the frequency 
errors by the mode method. In this case, K again approaches a constant value as the scale factor 
increases, but the particular constant value that it takes does depend on the SNR. As the SNR 
increases, however, the value of K eventually converges to a unique constant value below unity. 


Further explanation for this is given in the discussion below. 
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Fig. 7. (Color online) Top, closed symbols: The factor K =C“ / C , calculated directly from the master photon 
distribution converges to a constant value as the scale-factor enhancement S, increases. Bottom, open 
symbols: The factor K calculated from the frequency errors obtained from the mode method. The initial photon 
numbers (at S, =1, 7 =210K ) are 103, 10* and 10° (SNR values of 32, 100, and 1000), and decreases as S, 
increases. In the high scale-factor limit, K again becomes constant, but its value depends on SNR. 


Fig. 8(a) shows the peak frequency error, 6@ ),, , computed using the mean and mode methods 
discussed above, for the empty cavity (S, =1) and dispersive cavities (S, =344). The thick solid 
curve for the mean method is calculated from the RHS of Eq. (22), i.e., using the central limit 
theorem, and agrees with the simulation data. The peak frequency error for the two methods is 
identical when SNR=1. Then, as the SNR increases, the error for the mean method decreases 
rapidly and monotonically, whereas the error for the mode method undergoes an initial increase 
before eventually decreasing. This initial increase occurs because there are not enough photons to 
produce more than one photon inside any given frequency interval. Subsequently, both methods 
reach a plateau, where further increases in SNR are of less benefit. Note that the peak frequency 
error for the mode method is considerably higher than for the mean method. This is not a general 
result, but occurs in this case because the finesse is low. At higher values of the finesse, we find 
comparable errors for the two methods. This is expected because as the finesse increases the mean 


and mode of the distribution approach one another. 
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Fig. 8. (Color online) (a) QNL frequency error, 6@ 


‘ONL ? 
(bottom), calculated using the mean and mode methods. The abscissas are shifted with respect to one another 


at constant scale factor, S, =1 (top) and S, =344 


because for the same initial photon number, SNR is lower at the higher scale factor. Both methods eventually 
reach a plateau where further increases in SNR have little added benefit. (b) QNL enhancement, ¢2” for 
S, = 344 . The dashed line indicates the result obtained by setting K =1,i.e., the approximate result of Eq. (6) 
. The thick solid curves for the mean method are calculated from the RHS of Eqs. (21) and (22), i.e., using the 
central limit theorem. 


Fig. 8(b) shows the enhancement, ¢2%”, computed from the errors in Fig. 8(a) using both 


methods, at a fixed value of the scale factor (S, =344), as SNR is varied. We have assumed, in 
computing ¢2"”, that the method used for finding the transmission peak is the same for the 
dispersive and empty cavities. The dashed line in the figure represents the QNL enhancement 
predicted by Eq. (6), which falls between the upper and lower limits obtained from the mean and 
mode methods, respectively. Whereas the mean method yields results that are independent of SNR, 
the mode method produces enhancements that depend on SNR, and this dependency varies with 
the scale factor. However, at sufficiently high SNR, we find that the enhancement asymptotes to a 
finite (nonzero) value. This convergence to a nonzero value is a result of not being in the linear- 
dispersion regime. Consider that in the linear regime the mode broadens to infinite width and 
assumes a flat-top shape as the critical anomalous dispersion is approached. For a flat-top shape 


the error determined from the mode method would be finite even at infinite SNR. On the other 


hand, for the empty cavity (S, =1) the error tends towards zero as the SNR becomes infinite 


29 


because the mode shape is not flat-topped. Therefore, ¢°”” would asymptote to zero at high SNR. 


In contrast, we find that €2” converges to a nonzero value, which reflects the fact that the linear- 
dispersion approximation is violated at high scale factor. Although not shown in the figure, we 
have confirmed that the SNR required for this convergence increases with cavity finesse. 

The convergence to a finite QNL sensitivity enhancement is reiterated in Fig. 9(a), where the 
QNL enhancement £2”, computed using the mode method, is plotted against the scale-factor 
enhancement, S,. The mode method is not independent of SNR, rather it produces curves of ¢2"” 
vs. S, whose slope is lower for higher initial SVR values. Nevertheless, as the SNR increases, the 
slope eventually converges to the asymptotic value, corresponding to the flat region in Fig. 8(b). 
Fig. 9(b) shows the QNL sensitivity enhancement computed using both methods. Again, the thick 
solid curve is calculated from the RHS of Eq. (22), for the largest initial photon number, and agrees 
well with the simulation data obtained using the mean method, independent of the SNR. For the 
mode method, we only include the results for the highest initial photon number (N =10°) , where 
the slope no longer changes with SNR. 

Also plotted in Fig. 9(b) is the QNL enhancement calculated from Eq. (6), i.e., by setting K =1 
in Eq. (22), which again falls between the limits established by the two methods. Note at this value 
of K, the scale factor where the QNL enhancement crosses unity, i.e., where 2” >1, occurs at 
S, =6.4. On the other hand, this same crossing point, determined from the experiment described 
in section 6, occurs at S, =25 (see Fig. 4(b) and Fig. 6(b)). This difference is expected because, 
as pointed out in section 3, the scale factor in transmission is generally different from that in 
reflection owing to the differing amounts of reshaping that arise from the factor F. In addition, 


there is a slight difference in cavity finesse between the two cases due to the addition of the second 


output coupler. 
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Fig. 9. (Color online) (a) QNL enhancement, ¢°”, obtained from the mode method vs. scale-factor 


enhancement. At high scale factor, the curves are all linear. As the SNR increases the slope K decreases, 
eventually reaching the asymptotic value at high SNR. The initial (S, =1) photon numbers were N = 107, 103, 
10+, and 10° corresponding to SNR values between 10 and 1000. (b) QNL enhancement obtained from the mean 
(top, circles) and mode (bottom, squares) methods. For the mean method, different SNR values are represented 
by colored data points of varying color. For the mode method, only the asymptotic value is shown (initial 
photon number was N = 10°). The solid curve (top) is calculated from the RHS of Eq. (22). The dashed line is 
the result obtained from Eq. (6). 


Note, in Fig. 9(b), for both methods the enhancement initially decreases below unity, before 
eventually increasing proportionally with S,. This is because, as discussed in the main text, W, 
and M, approach finite values at the critical temperature. Therefore, in this high scale-factor limit 
(S, >> 1), the mode method produces QNL enhancements that are also described by Eq. (22), but 
with a constant value of K that differs from that of the mean method. The values of K calculated 
from the mode method for different initial SVR values are presented in Fig. 7, where they are seen 
to reach constant values at high S,. Moreover, while the particular constant value that K reaches 
depends on the initial SNR, a unique value is asymptotically approached at high SNR. Note that 
this procedure is lacking in comparison with application of the central limit theorem, which allows 
the frequency errors and QNL enhancement to be determined apriori by computing the constants 
C and K directly from the standard deviation and FWHM of the master photon distribution. 
Nevertheless, that the ratio K takes a unique constant value in the limits of high scale factor and 


SNR, signifies that the mode method can be used to readily establish a lower bound for the QNL 
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enhancement. In particular, once the value of K is established, knowledge of S,, M,, and W, is 
all that is required to characterize this lower bound. Finally, we note that in many experiments, 
e.g., in passive gyroscopes, the data range may be reduced far below one FSR. For such cases, the 
curves in Fig. 9(b) would tend to converge to the mode result, and the K =1 curve would fall 


outside the two bounds, revealing the inexactness of the approximation in Eq. (6). 


8. SUMMARY AND CONCLUSION 


Two distinct approaches were taken to investigate the possibility of obtaining a QNL 
sensitivity enhancement in the measurement of the mode frequencies of a passive cavity containing 
an anomalous dispersion medium. First, the scale-factor enhancement, normalized mode width, 
and signal-to-noise were simultaneously measured as the intracavity atomic medium is 
temperature tuned through the critical anomalous dispersion. Scale-factor enhancements larger 
than two orders of magnitude were obtained, translating to an order of magnitude enhancement in 
the predicted QNL measurement precision. A semi-empirical absorption model valid for low 
intracavity intensities was used to extract the temperatures and scale-factor enhancements. The 
accuracy of this model could be further improved by reducing systematic errors, i.e., through better 
design of the vapor cell oven and ensuring the proper empirical relationship between number 
density and temperature. Nevertheless, while these errors can shift the extracted temperatures (and, 
in turn, the critical temperature), they have negligible effect on the scale-factor enhancements. This 
is because, as we pointed out previously [6], the theoretical scale factor vs. temperature curve 
follows a ~1/T phenomenological relationship, regardless of whether the full treatment or the 
linear-dispersion approximation is used, or whether adjustments in the critical temperature occur 


as a result of correcting these systematic errors. Our model is, therefore, sufficiently robust to 
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extract scale-factor enhancements by fitting, even in the presence of systematic temperature errors 
(~1K). Moreover, the normalized mode width and signal-to-noise were determined directly from 
the cavity spectra by a model-independent process. Therefore, the enhancement in QNL precision, 
which is inferred from these quantities, is not significantly changed by these systematic errors. 

On the other hand, the assumption that we used to calculate the QNL error, tacitly ignored the 
effect of the change in mode shape, which is dramatic near the critical anomalous dispersion 
condition. To account for this, a second approach was taken that bypasses this assumption by 
computing the frequency errors directly through numerical Monte-Carlo simulations using the 
theoretical cavity spectrum as a master photon probability distribution function. The results of this 
numerical procedure were found to depend on the particular averaging method involved in finding 
the distribution peak. Nonetheless, we confirmed that the same product of S, W, and M can be used 
to deduce the QNL enhancement, after multiplication by a constant factor K whose value is found 
to fall within a set of bounds (see Fig. 7). 

Both approaches demonstrate the possibility of achieving an enhancement in measurement 
sensitivity in the QNL. To confirm the increase in precision, either in a passive ring gyroscope or 
linear configuration, the cavity would have to be held near QNL conditions, i.e., in a temperature- 
controlled vibration-free enclosure. This will be more difficult to achieve than for the equivalent 
empty cavity, in particular due to the temperature dependence of the atomic absorption [6]. It might 
be possible, however, to use a coupled-cavity, or other all-photonic design, rather than an atom- 
cavity configuration to reduce the temperature dependency. Such a confirmation would be an 
important proof-of-concept for the development of fast-light cavities and gyroscopes as a 


metrology tool. 
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